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■ A Weakly Interacting Massive Particle (WIMP) provides an attractive dark matter can- 

didate, and should be within reach of the next generation of high-energy colliders. We 
^ I consider the process of direct WIMP pair-production, accompanied by an initial-state 

5^ , radiation photon, in electron-positron collisions at the proposed International Linear 

Collider (ILC). We present a parametrization of the differential cross section for this 
process which conveniently separates the model-independent information provided by 
cosmology from the model-dependent inputs from particle physics. As an application, 
we consider two simple models, one super symmetric, and another of the "universal 
extra dimensions" (UED) type. The discovery reach of the ILC and the expected pre- 
cision of parameter measurements are studied in each model. In addition, for each of 
the two examples, we also investigate the ability of the ILC to distinguish between the 
two models through a shape-discrimination analysis of the photon energy spectrum. 
We show that with sufficient beam polarization the alternative model interpretation 
can be ruled out in a large part of the relevant parameter space. 



1 Introduction 



The most convincing and direct evidence for new physics beyond the Standard Model (BSM) 
is the existence of dark matter. The evidence for non-baryonic dark matter comes from a 
variety of observations on a wide range of astrophysical scales: rotation curves of galax- 
ies, cosmic microwave background, gravitational lensing, mass to luminosity ratios from the 
motion of galaxies in clusters of galaxies, power spectrum and large scale structure, and so 
on. No particle in the Standard Model (SM) can account for the observed dark matter. 
Interestingly, many extensions of the SM at the electroweak scale possess good dark mat- 
ter candidates. As the collider experiments begin directly probing this energy scale, it is 
important to understand what they can potentially teach us about dark matter. 

If one assumes that the dark matter particle was in thermal equilibrium with the rest of 
the cosmic fluid in the early universe, and then decoupled once the temperature dropped, 
its present abundance can be predicted. In a generic case, the abundance is a function of 
two continuous parameters: the dark matter particle mass and its total annihilation 
cross section into SM states at low velocity, (Tan- The current cosmological and astrophysical 
observations provide a measurement of the dark matter abundance, and indicate that the 
preferred values of (7an are of order 1 pb. Interestingly, this also happens to be the typical cross 
section for processes mediated by the SM weak interactions or other weak-scale particles. 
This surprising coincidence provides an additional hint that the dark matter particle may 
indeed by a part of the weak-scale theory. The dark matter particles of this class are called 
Weakly-Interacting Massive Particles, or WIMPs. 

There are two ways in which WIMPs can be produced in high-energy colliders. First, 
they can be directly pair-produced. (Since WIMPs are stable due to some new conserved 
quantum number, they cannot be singly produced in the collision of light SM particles.) 
Second, they can appear as decay products of other BSM particles. The second channel 
often has larger rates and smaller backgrounds, especially at the LHC where the production 
of strongly-interacting BSM states is enhanced. However, it is important to emphasize that 
this statement is very model-dependent and relies on the presence of other, heavier particles 
in the spectrum, which is a strong additional requirement independent of the WIMP dark 
matter hypothesis. In addition, the interpretation of the data in this channel is challenging. 
For example, the traditional methods for mass determination can provide a measurement of 
the WIMP mass only if long cascade decay chains are available and can be identified [1]. The 
measurement of the WIMP couphngs to the SM states is also a very challenging task and 
even in best-case scenarios one can only determine certain combinations of couplings and 
mixing angles [2] . Moreover, the WIMP may be only one of several weakly-interacting new 
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particles present in the model, and it may be quite challenging to disentangle the features 
of this sector in cascade decays. For these reasons, direct production of WIMPs, while more 
challenging experimentally, may offer unique opportunities to study the WIMP in the lab. 

The absence of a robust handle to distinguish events with direct WIMP production from 
the large SM background makes it virtually impossible to study this channel at the LHC 
[3, 4]. The relevant signatures (mono-jet plus missing energy and single photon plus missing 
energy) resemble direct graviton production in ADD models with Large Extra Dimensions 
[5], which is known to be among the experimentally most challenging searches at the LHC. 
Both channels are subject to a large irreducible background from Z + jets/photons, with 
Z decaying invisibly. In addition, there will be a large pure QCD background, where the 
missing energy arises as a detector effect. To make matters worse, the signal cross-sections in 
ADD are typically larger than the corresponding cross-sections for direct WIMP production. 
As a result of all these factors, the LHC is not expected to observe WIMPs directly, but 
only in the cascade decays of other, heavier particles. In this paper, therefore, we will 
consider the direct WIMP production at the proposed International Linear Collider (ILC), 
an e+e^ collider with = 500 GeV. Most previous studies have analyzed this process 
in the framework of supersymmctry (SUSY), and correspondingly refer to it as "radiative 
neutralino production" [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]. The 
related SUSY processes of direct production of sneutrinos or heavier neutralinos decaying 
invisibly, have also been considered [24, 25, 26, 27], leading to the so called "virtual LSP" 
scenarios. However, the connection between particle physics and cosmology provided by the 
WIMP hypothesis is very general, and is not restricted to supersymmctry. In particular, 
it may be readily applied in other model frameworks which have become popular in recent 
years, e.g. Universal Extra Dimensions (UED) [28, 29, 30, 31], Warped Extra Dimensions 
[32, 33, 34] and Little Higgs models with T-parity (LHT) [35, 36]. Given the multitude of 
existing models which can accommodate a successful WIMP, it is important to address the 
following two issues: 

1. When predicting the size of any potential signals of direct WIMP production at col- 
liders, how much of the discussion is model-independent, i.e. is relying only on the 
generic WIMP properties required by cosmology. 

2. Once a signal of direct^ WIMP production is observed, how can one discriminate among 
the various new physics model alternatives, e.g. SUSY, UED, LHT etc. 

^Model discrimination from indirect WIMP production at the ILC was previously discussed in [37, 38, 
39, 40]. 
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The first question was addressed in Ref. [3], which developed the model-independent 
formalism for predicting the rate for direct WIMP production at colliders. Using detailed 
balancing, the production cross-section at colliders was related to the leading term (Jan in the 
velocity expansion of the WIMP annihilation cross-section (the procedure will be reviewed 
below in Section 2). It was shown that the collider signals can be uniquely predicted in terms 
of only two continuous parameters: the WIMP mass and the WIMP annihilation fraction 
Kg into the e'^e~ final state (see the exact definition of Kg in eq. (7) below). Interestingly, 
it turned out that the ILC has sensitivity to direct WIMP pair-production in an interesting 
parameter range for {M^,Ke) [3, 41, 42, 43, 44], complementing more traditional methods 
for direct or indirect dark matter detection [43, 45]. 

While the formahsm of Ref. [3] is attractive due to its model-independence, it has one 
significant drawback: in order to use the connection to cosmology, the WIMP's must be 
produced at low velocity, i.e. near threshold. On the other hand, the center-of-mass energy 
of the ILC may very well be far above threshold. In order to bring the WIMP system back 
near threshold, a hard cut on the energy E^^ of the associated photon was required [3], which 
lead to a significant reduction in the observable signal. In principle, the photon energy cut 
may be relaxed [43], which results in larger signal rates, but then the connection to cosmology 
is lost and the accuracy of the theoretical prediction goes awry. 

In this paper, we tackle this problem by presenting a slightly modified framework to 
analyze this process. The new framework, outlined in Section 2, requires minimal additional 
inputs from the microscopic model of WIMPs, compared to the model-independent formalism 
of Ref. [3]. The advantage is that it allows to significantly relax the cut on the photon 
energy, extending the ILC reach. We also include the possibility of additional background 
suppression using polarized beams, further improving the reach [19, 20]. In Section 3 we 
show that once the observation is made, the ILC experiments can provide a measurement of 
the WIMP mass and the cross section of the process e+e~ — > XX- The latter measurement 
is especially interesting, since if this cross section (at threshold) turns out to be larger than 
that indicated by the WIMP relic abundance, the x particle will either be excluded as the 
sole dark matter component, or will require non-standard cosmology [46]. 

The second goal of this paper is to discuss the second question posed above: Once a 
signal of direct WIMP production is seen, how can one determine the WIMP properties 
and correspondingly the particular type of new physics model? This question has not been 
adequately addressed in the existing literature, since previous works on this channel have only 
considered one new physics model at a time, and did not attempt any model discrimination 
studies. In our case here, once our prediction of the observable photon spectrum becomes 
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model-dependent, a measurement of this spectrum can be readily used to distinguish between 
different models of WIMPs and their interactions. The results from such an analysis will be 
presented in Section 4. 



2 Direct WIMP Production at the ILC 

In this section, we will define the framework to describe a simple observable signature of the 
direct WIMP production in electron-positron collisions, and estimate the expected reach of 
the ILC in this framework. 

2.1 The Signature 

We consider the direct WIMP pair-production process 

e+e" XX- (1) 

Here, the WIMP x is an electrically neutral, color-singlet particle, with a mass in the 1 
GeV-few hundred GeV range. At this point, we do not need to specify any other quantum 
numbers of x (such as its spin), or the fundamental theory of which x is a part. On its own, 
the process (1) docs not leave an observable signature in the detector. To get an observable 
signature, we consider the closely related process 

e+e~ XX7, (2) 

in which the WIMP pair can be observed as missing momentum recoiling against the detected 
photon. The cross section is dominated by the photons that are soft or coUinear with the 
beam line: 

EySine <^ ^/s , (3) 

where is the photon energy, ^/s is the center-of-mass energy of the collider and 6 is the 
angle between the photon and the incoming electron. In this regime the differential cross 
section can be factorized according to the Weitszecker- Williams formula: 

where = a{e'^e'' XXT); ^2 = cr{e'^e~ — > xx)j ^ — '^^-y/V^^ ^ — •^(1 — x), and 

^(,,,).i^i±(i^ 1 . (5) 

TT X sm 9 
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Of course, very soft or coUinear photons remain undetected in a realistic detector; however, it 
was shown in Ref. [3] that the factorization (4) works rather well under realistic experimental 
conditions. 

In general, the function (72(5) depends on the microscopic model of the WIMP and its 
interactions with electrons. However, in the region close to production threshold, s — 4M^ <^ 
s, this function can be approximated as 

a2{s)^aoil-^\ , (6) 

where p is the angular momentum of the lowest partial wave contributing to this process. (In 
particular, p = for s-wave scattering and p = 1 for p-wave.) This is the approximation that 
was used in Ref. [3]. Its advantages are the independence from the underlying microscopic 
model, and the fact that the quantity ctq can be easily related to the relic abundance of the 
X particle. The key relation is [3] 

ao = 2'(P-'^{2S^ + iyK,a,^, (7) 

where is the spin of the x particle, (Tan is the coefficient of the leading term in the small- 
velocity expansion of the total x pair-annihilation cross section, and Kg is the fraction of 
X pair-annihilation events at low velocity which result in electron-positron pairs. At this 
point we also assume that x = X: that the WIMP x is identical to its antiparticle. This 
assumption will be true in our two examples discussed explicitly below, but is not really 
necessary in general - in the cases where x 7^ x, there is simply an additional factor of 2 in 
the right-hand side of (7). 

The quantity is directly related to the present relic abundance of the x particle. 
Assuming that x makes up all of the observed dark matter, and using the observed value for 
the dark matter density fldmh'^ = 0.1143 ±0.0034 [47], yields the determination of o"an shown 
in Fig. 1. The quantity cannot be obtained from cosmology, and needs to be calculated 
once the model of particle physics is specified. Note however that in any model Kg < 1, so 
that eq. (7) provides a model-independent upper bound on ctq: 

<J0 < 22(*'-2) (25^ + 1) Van , (8) 

which is valid as long as x is the dominant dark matter component. If the value of ctq 
measured in collider experiments turns out to violate the bound (8), there can be only 
two, equally exciting, explanations: that in addition to the x particle, there must exist 
yet another independent WIMP candidate, or that there was some sort of non-standard 
cosmological evolution in the early universe. 
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Figure 1: Values of the quantity a an allowed at the Sa level as a function of the WIMP mass 
M^, and for different values of the WIMP spin S-^. The lower (upper) band is for models 
where s-wave (p-wave) annihilation dominates. We use the constraint Jldm^^ = 0.1143 ± 
0.0034, which results from the combination of data from WMAP-5, Type la supernovae and 
baryon acoustic oscillations [47]. 



The main disadvantage of the approximation (6) is its hmited range of vahdity: if one 
wishes to use it to analyze data, one needs to impose a lower bound on the photon energy to 
restrict the analysis to the near-threshold region, thus cutting out most of the signal. In this 
paper, we will pursue a different approach: we will calculate (72(5), for all values of s, in two 
benchmark models, and describe the ILC sensitivity for these models. As we will see, this 
gives a substantially better sensitivity than the model- independent search of [3] . To preserve 
the simple connection to relic abundance, we choose to parametrize the cross section as 

r2 \ p+1/2 



0-2(s) = (^of{s) 1 - 



4:Mi 



(9) 



where f{s) is a model-dependent dimensionless shape function obeying the normalization 
condition 

/(4M^) = 1 . (10) 

This form reduces to eq. (6) in the near-threshold region, so the coefficient ctq is related to 
the relic abundance as described above. 

Note that effects such as beamstrahlung [48] and the possibility of multiple ISR emission 
are not included in our approximation scheme. In the future, it would be interesting to 
incorporate these effects in the analysis. 



6 



X 




e+ X e 

(a) (b) 



X 



X 



-X 



■X 



(d) 



Figure 2: Feynman diagrams for the WIMP pair-production in the SUSY benchmark model 
(a,b) and the UED benchmark model (c,d). 

2.2 Benchmark Models 

Our first benchmark model is supersymmetric [49], with a pure Bino LSP pair-produced in 
electron-positron collisions via t-channel exchange of the right-handed selectron ept, see the 
diagrams in Fig. 2 (a,b). This simple toy model approximates the constrained MSSM in the 
"bulk" dark matter region, up to subleading contributions due to non-bino admixtures in the 
neutralino and the heavier cl exchange. In this model, p = 1, since s-wave annihilation at 
threshold is forbidden by CP invariance (up to terms suppressed by m^/M^ ~ lO^^*', which 
we ignore). We obtain 
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Note that Q < y, z < 1, and that Ymiy^i f {y , z) = 1, as required by eq. (10). This is the 
cross section for unpolarized beams; since only right-handed electrons couple to en, the 
corresponding cross section for polarized beams can be found by simply multiplying by an 
overall factor of 

(l + ^e-)(l 

where 



(13) 
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are the electron and positron beam polarizations, respectively. 

The second benchmark model is a spin-1 WIMP, pair- produced in electron-positron col- 
hsions via i-channel exchange of the spin- 1/2 "heavy electron" Er, see the diagrams in 
Fig. 2(c,d). The heavy electron is a Dirac fermion, and is assumed to be the partner of 
the right-handed electron cr, this assumption fixes the helicity structure of the vertices 
in Fig. 2(c,d). This toy model reproduces the leading contribution to this process in the 
minimal UED model [29, 50, 51].^ In this case, s annihilation is allowed, and p — 0. We 
obtain 



0-0 



f{y,z) 



1287rM|^ (1 + ^)2' 



where 



128z^^/T^{y{l- z) + 2z){y{l- zy + 4:z 
4:z^l-y(y{l -z) + 2z) (y{l - zf{Sz^ + 3) + Az{Sz^ - Az + 3)) 
- - zf{Sz^ + 3) + Ayz{8z^ - 6z + 3) + 8z^{Az'^ + 1)) 

y{l - zf + 4.) log ^[^-^l + ^f -fZZ;] , (15) 
^ ' ^ ^z/(l-^)+2^(l + VT^)J' ^ ' 

AMI Ml ^ ^ 

^^-T' '^Jit- '''' 

Note that we will continue to use the "generic" WIMP notation, x, in the UED case; the x 
in this case is actually the spin-1 massive partner of the hypercharge gauge boson, commonly 
known as [53]. 

Of course, it is straightforward to apply our formalism to other models, and we expect 
qualitatively similar results, as long as WIMPs do have an unsuppressed coupling to electrons. 
In all such cases the starting point of the analysis would still be the parametrization (9), 
from where the current discussion can be generalized in two different ways. First, one may 
consider a different region of parameter space of the same model (SUSY or UED), where 
there can be additional relevant diagrams beyond those shown in Fig. 2. For example, the 
addition of a i-channel cl exchange would simply modify the functional form of (11) and 
introduce an extra parameter (the mass of e^,) . Alternatively, one may consider a different 
model altogether, by deriving the corresponding function /(s) relevant for that case. From 



^The unpolarizcd cross section is identical in the minimal Littlest Higgs model with T parity [52, 36], 
but the helicity of the vertices is the opposite, since only the left-handed electron has a partner in that 
model. Note also that the value of Kg in the LHT model is small, of order 10~*, due to a group theory factor 
suppression of the vertices in Fig. 2(c,d). 
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that point on, the analysis would still follow the steps outlined below for our two model 
examples (SUSY and UED). 

2.3 ILC Reach 

The above formalism provides a prediction for the differential rate of 7+missing energy 
events at an e^e~ collider in terms of two continuous parameters, M-^ and ctq, once a specific 
model is chosen to fix the function /(s). Experiments at an e^e~ collider can search for this 
signature as an excess over the SM background. We analyzed the reach of a 500 GeV ILC, 
assuming the following kinematic acceptance cuts: 

sin ^ > 0.1, pt,7 = E^sme>2> GeV. (17) 

The pr,7 cut is imposed to reject backgrounds such as e'^e" — > e^e"^ in the region where the 
and e~ are too forward to be detected. This cut corresponds to the improved BeamCal 
acceptance of 6.67 mrad for 500 GeV hnear collider [54]. With these cuts, the SM background 
is strongly dominated by the process e'^e~ ^i^J, and we ignore all other backgrounds in 
this study. We simulated this background using CalcHEP [55] and CompHEP [56]. To quantify 
the reach, we compare the differential cross sections da /dE^ for signal+BG and pure EG. To 
account for finite detector resolution, we smear the energy of the photon spectrum for both 
signal and background according to 6E/E = 14.4%/v^gcV © 0.5% [41]. We bin the events 
in 5 GeV bins in the photon energy, and compute the event numbers in each bin, {iVsig+Bc) 
and {N^q}, for a given integrated luminosity £int. We then use the standard technique 
to quantify the probability of measuring the event numbers predicted by the BG-only model 
if the true model is signal+BG. We define 

where 5sys is the fractional systematic uncertainty in the background prediction. Note that 
we conservatively assume that this systematic error is uncorrelated between bins. 

The reach of the search in the benchmark models defined above is presented in Fig. 3. 
We assumed an integrated luminosity of C\nt = 500 fb^^, and a systematic error on the 
background prediction 5sys = 0.3%. The 3(T evidence and 5cr discovery contours correspond 
to = 9 and 25, respectively. The contours labeled "polarized" correspond to electron and 
positron beam polarizations of Pg- = 0.8 and Pe+ — 0-6, correspondingly. Note that these 
polarizations strongly suppress the dominant contribution to the background from i-channel 
W exchange, which only appears for left-handed electrons, while at the same time enhancing 
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Figure 3: The reach of a 500 GeV electron-positron colhder with an integrated luminosity of 
^int = 500 fb~^ for discovery of our two benchmark models: (a) SUSY and (b) UED. Blue 
(red) lines correspond to unpolarized (polarized with P^- = 0.8, Pe+ = 0.6j beams. Solid 
(dashed) lines indicate 5a (3a ) significance. Also shown are Ke contours in the (a) (M^, Mg^) 
or (b) {M^,Meji) parameter plane. 



the signal in both of our benchmark models, where only the right-handed electrons couple to 
WIMPs. The reach is presented in terms of the WIMP mass and the mass of the t-channel 
particle which dominates the production process. One should keep in mind that the t-channel 
particle can also be directly produced at the ILC, provided its mass is within the kinematic 
reach of the collider. In this sense, one should note that the 5a reach with unpolarized beams 
in the SUSY case of Fig. 3(a) only includes regions where the right-handed slepton e/j would 
be directly produced as well. Consequently, in order to access the remaining parameter space 
where the radiative neutralino production is the only available signal, one must rely on beam 
polarization. In Fig. 3 we also show contours of constant Kg, defined in eq. (7). For low 
WIMP masses, the reach of the ILC is impressive: with polarized beams, the WIMPs below 
100 GeV can be discovered even if the e~^e~ channel only contributes 0.1% to the total WIMP 
pair-annihilation cross section in the early universe. For larger WIMP masses, the sensitivity 
decreases due to the kinematic suppression in WIMP production. This suppression can of 
course be compensated by a higher center-of-mass energy. 
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3 Parameter Measurement 



Once the presence of non-SM contribution in the 7+missing energy channel is estabhshed, 
one can attempt to fit this contribution with BSM model predictions. In each model, the 
prediction depends on a number of unknown parameters; the fit to data can be used to 
constrain these parameters. In this section, we will discuss the accuracy of this procedure, 
using our benchmark models as examples. (Of course, it is possible that in a given model 
a good fit cannot be obtained for any sensible parameter values; in this case, the model 
would be ruled out as the interpretation of the data. We will discuss examples of such model 
discrimination in the next section.) 

In the first example, we assumed that the data is given by the prediction of the SUSY 
benchmark model (plus SM background), with 

= 100 GeV, Mg^ = 300 GeV. (19) 

For reference, these values correspond to Ke — 0.023. We then performed a scan of the 
model parameter space, varying M-^ between 50 and 250 GeV, and Mg^ between 100 and 
1000 GeV. (We imposed Mg^ > M^.) For each point in the scan, we computed the probability 
of observing the data given the model parameters at that point. To quantify this probability, 
we used the same technique as in the reach analysis above. The results are presented 
in Fig. 4(a). With unpolarized beams, only a very crude mass determination is possible: 
for example, any WIMP mass up to 220 GeV is allowed at a 2-sigma level. With polarized 
beams, however, a rather impressive accuracy of about ±20 GeV in both and Mg^, at a 
2-sigma level, can be achieved. 

In Fig. 4(a) we also show contours of the quantity Kg, which is uniquely fixed by eq. (7), 
once the values of and Mg^ are specified. Therefore, a measurement of the parame- 
ters (M^, Mg^) can be immediately reinterpreted as a n^. measurement. In particular, from 
Fig. 4(a) we see that with unpolarized beams, the ILC will only set an upper bound on 
of about 0.22 at the 3(t level, while using the polarized beam option, the ACg determination 
is much more precise: 0.01 < Ke < 0.04 at 3(7. 

In the second example, we repeated the same exercise for the UED benchmark model. 
The results, presented in Fig. 4(b), are qualitatively similar. The main difference in this case 
is quantitatively better sensitivity, mainly due to a higher signal cross section for the chosen 
model point (M^ = 100 GeV, M^^ = 300 GeV, which in this case corresponds to k^. — 0.53). 
As a result, a fairly accurate measurement of the masses is possible even with unpolarized 
beams: the accuracy on the WIMP mass is of order ±10 GeV (±20 GeV) at the 2-sigma 
(3-sigma) level. With polarized beams, a WIMP mass determination at the level of ±2 GeV 
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Figure 4: Expected precision of the model parameter measurement in the two benchmark 
scenarios: (a) SUSY and (b) UED. The notation and labeUing are the same as in Fig. 3. 



becomes possible. The corresponding 2-sigma determination of is now 0.2 < Kg < 1.3 
(0.46 < < 0.6) with unpolarized (polarized) beams. 

Note that the capabihty of the ILC to measure the WIMP mass in this process has also 
been recently analyzed by Bartels and List [44]. Their results for the spin-1 case (the only 
one analyzed in their study) are quahtatively in agreement with ours. Precise agreement is 
not expected, since the two studies used a different choice of the benchmark point, as well as 
a slightly different statistical analysis procedure - for example, Ref. [44] fits the WIMP mass 

assuming that Kg is already known, whereas we perform a 2-parameter fit. Our results are 
also in qualitative agreement with those of Refs. [42, 43], which used the mo del- independent 
framework of [3] to perform a 2-parameter fit for (M^, Kg). 



4 Model Discrimination 

Since the shape of the photon spectrum in the events with WIMP production is model- 
dependent (as characterized by the function f{s)), measuring this shape can provide inter- 
esting information about the underlying model. While not an unambiguous spin determina- 
tion, this measurement can nevertheless distinguish scenarios in which WlMPs have different 
spins, as long as the WIMP interactions are specified in each model. As always, it is impor- 
tant to keep in mind that each of the models has free parameters, and a model can be "ruled 
out" only if the point in its full parameter space providing the best fit to data still does not 
provide an acceptable fit. (This point has been recently emphasized in Refs. [57, 2, 58].) 
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Figure 5: Plots illustrating the ability of the ILC to discriminate between the two benchmark 
scenarios using the radiative WIMP production process. 

To illustrate the power of the ILC 7+missing energy measurements to distinguish between 
the two benchmark models, we performed two sample studies. In the first one, we assume 
that the data is given by the prediction of the SUSY benchmark model (plus SM background), 
with the same parameters as in eq. (19). We then attempt to fit the data with the predictions 
of the UED model, varying the parameters and Me^- For each parameter point, we 
determine the value, and then compute the conditional probability p that the UED 
model with these parameters is true, given the data. The results are shown in Fig. 5(a). 
With unpolarized beams, the model-discriminating power of this measurement is rather poor: 
any UED model with parameters above the upper blue (dashed) contour could be true with 
probability above 40%. The UED model needs to be in the high- mass region to match the 
observed overall event rate, but since this rate is rather low and the background is high, there 
is not much sensitivity in this case. The model-discriminating power is greatly enhanced by 
beam polarization, which strongly suppresses the background and (in both of our benchmark 
models) enhances the signal. With polarized beams, the UED interpretation can be safely 
ruled out: even the best-fit point in the UED space (which happens to be at = 140 GeV, 
Meji = 621 GeV) has only a 0.5% probability of describing the data. The photon spectra for 
the data (with polarized beams) and the best-fit UED point are shown in Fig. 6(a); it is clear 
from this figure that the quality of the UED fit to data is rather poor. One can conclude 
that the UED benchmark model as a whole is ruled out at a 99.5% CL by this data. 

We repeated this analysis assuming that the true model is UED, with 

= 140 GeV, Me^ = 621 GeV, (20) 
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Figure 6: Signal photon spectra in the two benchmark models, (a) Black histogram: SUSY 
model with the parameters in eq. (19); red histogram: UED model with parameters giving 
the best fit to the SUSY model, (b) Black histogram: UED model with the parameters in 
eq. (20); blue and red histograms: SUSY models with parameters as indicated by the two 
crosses in Fig. 5(a). In both cases, polarized electron and positron beams are assumed. The 
error bars on these ''signal" histograms are derived from the full signal+BG distributions, 
with the statistical (Cint = 500 fb"^) and systematic (0.3%) errors, added in quadrature. 

and attempting to fit it with the SUSY predictions. The results are shown in Fig. 5(b). 
In this case, the model-discriminating power of the measurement is significantly weaker, 
primarily due to the low signal cross section. In the unpolarized case, most of the SUSY 
parameter space is allowed; only the points with low and Mg^ are constrained. The 
situation is improved with the polarized beams, where most of the SUSY parameters space 
is ruled out at a confidence level of 99.9% or better. However, a significant part of the 
parameter space survives. This is clear from the photon spectra shown in Fig. 6: the "data" 
spectrum is well fit, within the error bars, by the predictions of sample SUSY points. The 
best-fit SUSY point, at = 170 GeV, M^;^ = 200 GeV, gives a photon spectrum which is 
essentially identical to the input UED "data" . Thus, convincing discrimination between the 
two models appears impossible in this case. However, this rather pessimistic conclusion is 
largely due to our choice of a relatively heavy UED study point. When we repeat the same 
analysis for the UED study point used in Fig. 4(b), namely = 100 GeV and Mej^ = 300 
GeV, we find that the SUSY interpretation is now completely ruled out. 

The two models selected for our study differed in the behavior of the cross sections in the 
non-relativistic limit (s-wave for UED vs. p-wave for SUSY). This leads to a large difference 
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in the normalization factors ctq in the two models, see Eq. (8). However, by varying the 
masses, the overall rates of the two models can be made to agree, and it is the photon spec- 
trum shape difference that is crucial for breaking that degeneracy. If instead we considered 
two models that both had, for example, s-wave annihilation, then the normalization ctq would 
be the same, and one would need to rely exclusively on the spectrum shape differences. As 
long as the spectra predicted by the two models are sufficiently different, model discrimi- 
nation should be possible. Thus, the proposed technique should be useful in discriminating 
between models with the same threshold behavior of the WIMP annihilation cross section. 

5 Conclusions 

The ILC may offer a unique window into the dark matter world, by providing an opportunity 
to directly produce dark matter particles and study their properties. In this paper we revisited 
what is perhaps the most conservative and least model-dependent collider signature of dark 
matter: direct pair production of WIMPs, tagged with a photon from initial state radiation 
[3]. Our main results are as follows: 

• We extended the formalism of Ref. [3] into the region away from the endpoint of the 
photon energy spectrum. The advantage of this new, although somewhat more model- 
dependent, approach is that it utilizes all of the available information from the single 
photon data, and not just near the endpoint. As a result, the discovery reach of the 
ILC is significantly improved, since one is going beyond a simple counting experiment, 
and using the measured shape of the spectrum as well. 

• In Section 2.3 we estimated the ILC discovery reach for two representative models: 
SUSY and "UED"-hke. In both cases we find that radiative WIMP production can be 
observable, as long as the WIMP mass is within the kinematic reach of the ILC and 
the WIMPs have non-negligible couplings to electrons. It is worth emphasizing the 
large difference in the reach with and without beam polarization - see Fig. 3. 

• As demonstrated in Sec. 3, the ILC can determine the masses of both the WIMP and 
the associated t-channel particle (in our two examples, a slepton or a KK-lepton). This 
measurement is particularly impressive if beam polarization is available - see Fig. 4. 

• Having observed a signal of radiative WIMP production, the ILC may also attempt a 
discrimination among different model explanations. The outcome of this analysis in 
general depends on the specific region of parameter space where the data happens to 
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lie, but as a rule, beam polarization significantly improves the chances for success - 
see Fig. 5. 
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